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We systematically apply Ginzburg-Landau theory to determine BCS pairing in a strongly-coupled 
uniform superfluid of three-flavor massless quarks in flavor equilibrium. We elucidate the phase dia- 
gram near the critical temperature in the space of the parameters characterizing the thermodynamic- 
potential terms of fourth order in the pairing gap. Within the color and flavor antisymmetric channel 
with zero total angular momentum, the phase diagram contains an isoscalar, color-antitriplet phase 
^~>) ■ and a color-flavor-locked phase, reached by a second order phase transition from the normal state, 

' as well as states reached by a flrst order phase transition. We complement the general Ginzburg- 

Landau approach by deriving the high-density asymptotic form of the Ginzburg-Landau free energy 
from the flnite temperature weak-coupling gap equation. The dynamically-screened, long-range 
) ^ color magnetic interactions are explicitly taken into account in solving the gap equation. We find 

. that in the limit of weak coupling, the isoscalar, color-antitriplet phase has higher free energy near 

QQ ' the transition temperature than the color-flavor locked phase. 

In view of the fact that deconfined quark matter must be color charge neutral, we incorporate 
the constraint of overall color neutrality into the general Ginzburg-Landau theory and the gap 
equation. This constraint yields a disparity in the chemical potential between colors and reduces 
^ , the size of the pairing gap, in the presence of the anisotropy of the order parameters in color space. 

In comparison with the case in which there are no chemical potential differences between colors and 
hence the superfluid generally has nonzero net color charge, we flnd that while the constraint of 
color neutrality has only negligible effects on the gap in the weak coupling regime, it appreciably 
destabilizes the isoscalar, color-antitriplet phase in the strong coupling regime without affecting the 
color-flavor-locked phase. 
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I. INTRODUCTION 

The possibility that degenerate relativistic quark matter becomes a color superconductor at low temperatures has 
, . • been considered for the past two decades. As the seminal papers noted, the quark-quark interaction in the 
r> [ color antitriplct channel is attractive and drives a Cooper pairing instability in the system, even in the limit of high 
1 densities where the Fermi energy of the quarks dominates over the one-gluon exchange interaction energy. Work in 
" " ' the intervening period has concentrated on the mean field approach, strictly valid only in the weak coupling limit 
(see, e.g., Refs. [|| and |^ for reviews). As the density is lowered, nonperturbative efi^ects arising from self-couplings 
of the gluon field prevail, finally leading to a confinement transition into hadronic matter. Strong coupling effects can 
modify the equilibrium order parameter from the weak coupling prediction, as effects of spin fluctuation exchange do 
in superfluid '^He The resultant change in the color-superconducting phase can also affect the breaking of chiral 
symmetry |^-|^. The properties of the color superconductivity such as the pairing gap and the critical temperature, 
Tc, have yet to be derived in the strong coupling regime in a rigorous way. 

In this paper we examine color superconductors by means of a general Ginzburg-Landau approach, which permits 
us to determine the coarse grained features of such systems at temperatures just below Tc, for arbitrary QCD coupling 
constant, g. This approach, pioneered by Anderson and Brinkman ||^ and Mermin and Stare pO| ] in the context of 
superfluid '^He, reveals the most energetically favorable phase just below Tc, in terms of the parameters characterizing 
the thermodynamic potential to fourth order in the pairing gap. The first application of the general Ginzburg-Landau 
theory was made by Bailin and Love ||^ to quark pairing with one flavor and total angular momentum J — 1. Here 
we shall consider more general Cooper pairing between u, d, and s quarks with J = 0. We complement this approach 
by deriving the parameters controlling the fourth order terms from the weak coupling gap equation. Previous work 
has systematically investigated such pairing at zero temperature; the present study provides a systematic 
elucidation of the equilibrium properties near Tc. 
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The zero-temperature pairing gap and hence Tc is predicted to be 10-100 MeV for baryon chemical potentials ~ 
1 GeV [p^ . This prediction relies on extrapolation from weak coupling into the low density, nonperturbative regime, 
by incorporating into the BCS gap equations effective interactions modeled after instanton-mediated interactions, in 
such a way as to reproduce constituent quark masses. Application of the Wilson renormalization group to analysis of 
the stability of a Fermi liquid against Cooper pairing |l^ suggests that four-fermion couplings, induced between two 
particles with zero total momentum by one-gluon exchange or by instantons, grow logarithmically with momentum 
as higher modes are successively integrated out closer to the Fermi surface; the scattering amplitudes eventually 
reach a singularity - a Landau pole - in a way dependent on the number of flavors involved in the pairing. For two 
flavors (u, d), both gluon and instanton-induced interactions play a role in opening an energy gap in the isoscalar, 
color- antitriplet channel with zero total angular momentum (J = 0). Pairing in this channel partially breaks baryon 
number symmetry and global color rotational invariance, but restores chiral symmetry. For three flavors (u, d, s), on 
the other hand, a color-flavor locked state |l^ arises in the J — Q channel mainly from the gluon-induced interactions 
when the strange quark mass is sufficiently small. This state is invariant under simultaneous exchange of color and 
flavor, but not under single exchange of color or flavor; chiral symmetry as well as baryon number symmetry and 
global color rotational invariance are broken. In the high density regime where the interactions are dominated by 
one-gluon exchange, the color magnetic (transverse) force, which is screened only dynamically by Landau damping 
of its mediators |l|,|l^ , is sufficiently long ranged to alter the dependence of the pairing gap on the QCD coupling 
constant g from the BCS result. This fact was first noted by Son ||l^ using a renormalization group method and an 
approximate solution to the relevant gap equation. 

As in prior papers, we focus on the equilibrium properties of an ultrarelativistic color superconductor that is 
homogeneous, in the sense that the real gluon field vanishes everywhere and the order parameter is everywhere 
continuous in magnitude and orientation; we denote such a state as superfluid quark matter. Such homogeneity is 
similar to that in superfluid '^He and superfluid neutron matter, as noted by Bailin and Love because in both cases 
breaking of the global U{1) gauge symmetry is accompanied by global symmetry breaking associated with the internal 
degrees of freedom. In superfluid quark matter, the possible order parameters are generally anisotropic in color space 
(see, e.g., Ref. |2)), a situation analogous to superfluid "^He in which, as seen experimentally, the anisotropy lies in 
spin space 

We first restrict ourselves to the case, normally assumed in earlier investigations, in which the chemical potentials of 
different color and fiavor quarks are equal. We obtain the thermodynamic potential difference between the superfluid 
and normal phases near T^, from the Ginzburg-Landau approach. The terms of second and fourth order in the pairing 
gap are constrained by invariances of the grand canonical Hamiltonian and by the structure of the order parameters, 
assumed here to be antisymmetric in color and flavor space. We then identify the degenerate sets of order parameters 
corresponding to local energy minima as isoscalar color-antitriplet and color-flavor locked states, and determine 
their condensation energies. We draw the resultant equilibrium phase diagram near Tc in the space of parameters 
characterizing the fourth order terms. We find that in the limit of weak coupling, the isoscalar, color-antitriplet phase 
is less favorable than the color-fiavor locked phase near Tc. 

In general, determination of the Ginzburg-Landau parameters requires inclusion of strong coupling effects. In the 
weak coupling limit, we anchor the general Ginzburg-Landau approach by deriving the parameters from the relevant 
weak coupling gap equations at finite temperature P,|l^,^, including the infrared structure of the gluon propagator 
as in Ref. [ pO[ . The latter behavior determines the g dependence of the pairing gap in the weak coupling limit. Debye 
screening of the electric gluons and Landau damping of the magnetic giuons, calculated with one-gluon exchange as 
modified by a normal medium, provide effective infrared cutoffs of the respective scattering amplitudes. 

Motivated by the fact that deconfined quark matter must be in an overall color singlet state, we also consider the 
equilibrium properties of superfluid quark matter under the condition that it is color-charge neutral. In normal quark 
matter, in which global color rotational invariance is ensured by QCD interactions, the constraint of color neutrality 
leads to equality of the chemical potentials for quarks of different colors (see, e.g., Ref. [^). Superfluid quark matter, 
on the other hand, generally has a preferred direction in color space. Such violation of color rotational invariance 
together with the requirement of overall color neutrality leads to differences in the chemical potential between colors, 
in a way dependent on the degree of anisotropy of the order parameters in color space. This disparity in chemical 
potentials in turn alters the pairing gap, a feature seen in the general Ginzburg-Landau analysis. Beyond the usual 
terms of second and fourth order in the pairing gap, color neutrality adds terms dependent on the chemical potential 
differences between colors. Such extra terms, which act in the isoscalar, color-antitriplet channel, do not remove the 
degeneracy of the order-parameter sets occurring for equal color chemical potentials and renormalize the coefficients 
of the fourth order terms so as to reduce the pairing gap. In order to estimate the influence of color neutrality on 
the gap in the weak coupling limit, we incorporate color chemical potential differences in the gap equation, and flnd 
that in the weak coupling limit, such differences do not significantly modify the gap. Although it does not affect 
the color-flavor-locked phase, the constraint of color neutrality can appreciably raise the free energy of the isoscalar, 
color-antitriplet phase in the strong coupling regime. 
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We base our calculations in the weak coupling limit on the field-theoretical description of an ordinary supercon- 
ductor as advanced by Nambu |^ and Eliashberg who reformulated the original BCS problem in terms of 
electron-phonon interactions. Such a formalism is useful when the interactions are nonlocal in time, as are the color 
magnetic interactions in the relativistic plasma. This tool was generalized to a relativistic regime by others in earlier 
investigations (e.g., Ref. |^). 

We focus in this paper on a system in flavor equilibrium, with no chemical potential differences between flavors (as 
long as the electrical charge of the quarks is neglected) . As a consequence, the energy of the superfluid is invariant 
under rotations of the equilibrium states in flavor space. The influence of the possible flavor chemical potential 
differences on the phase diagram, examined by others in earlier investigations 0,|2^, resolves this degeneracy. Such 
differences may also give rise to an inhomogeneous phase in which the pairing gap varies periodically in ordinary 
space. The conditions determining the chemical potential differences between flavors depend, however, sensitively on 
the extent to which the system is out of flavor equilibrium. 

We summarize the main conclusions of this paper. By applying general Ginzburg-Landau theory to the J — 0, 
color and flavor antisymmetric pairing states of superfluid quark matter, we elucidate the phase diagram near Tc in 
the space of the parameters controlling the fourth order terms in the pairing gap. We find that the phase diagram 
contains an isoscalar color-antitriplet phase and a color-flavor locked phase, both reached by a second order phase 
transition, as well as a region of Ginzburg-Landau parameters for which the transition to the superfluid phase is of 
first order. The detailed structure of the superfluid phase in this regime depends on the nature of the sixth order 
terms in the expansion. In the limit of weak coupling, the color-flavor locked phase is more favorable than the isoscalar 
color-antitriplet phase. We also expand the general Ginzburg-Landau approach and the Nambu-Eliashberg formalism 
to incorporate differences in the chemical potential between colors, required to preserve color-charge neutrality of the 
system. We find that the color neutrality constraint in the presence of anisotropy of the order parameters in color 
space tends to suppress the gap. 

Color superconductivity in quark matter becomes an astrophysically interesting problem if neutron star interiors 
are sufficiently dense that they contain quark matter cores p5| . Generally, a quark superfluid in a neutron star would 
not be electrically neutral since each quark has fractional electric charge; rather it would coexist with electrons (and 
muons) in such a way as to ensure electric neutrality in the system. Due to the dually charged nature of the quarks, 
macroscopic manifestations of both color and electromagnetic superconductivity such as Meissner effects, generation 
of London fields, and vortex formation are expected from magnetic fields and rotations as observed in these celestial 
objects. We will discuss these issues, which may be relevant to magnetic structure, cooling, and rotational evolution 
of the neutron stars |^,^ , elsewhere |^ . 

In Sec. 1^ we construct the generalized Ginzburg-Landau theory including color chemical potential differences, and 
apply it to color and flavor antisymmetric pairing with J — 0. The equilibrium phas e d iagrams for various values 
of the parameters controlling the expansion of the free energy are given here. Section |lll| is devoted to deriving the 
gap equations in the weak coupling limit, into which the color chemical potential differences are incorporated, and to 
calculations of the energy gap for the two types of optimal pairing states derived in Sec. In Sec. IV, we calculate 
the thermodynamic potential difference between the superfluid and normal phases near Tc from the structure of the 
energy gap obtained in Sec. Ul . Our conclusions are given in Sec. M. We use units in which ?i = c = fcs = 1. 



II. GENERAL GINZBURG-LANDAU APPROACH 



In this section, we apply the general Ginzburg-Landau approach to a quark superfluid that is in an overall color 
singlet state. We then derive the thermodynamic potentials in terms of the color chemical potential differences. We 
finally draw the equilibrium phase diagrams near Tc for the parameters characterizing the strong-coupling effects and 
the color-neutrality constraint. 

We consider uniform quark matter of temperature T and baryon chemical potential /i, with number of flavors 
Nf = 3 {u,d,s) and colors Nc — 3. We neglect quark masses and electrical charges as well as chiral condensates, 
and do not take into account at this stage either the color or electromagnetic gauge field, except for virtual gluons 
mediating interactions between the quarks. We also assume that the system has zero net color charge and satisfies 
the condition for flavor equilibrium, 

Mau = Mad = A^as = Ma , (1) 

where iiai is the chemical potential for the quarks of color a and flavor i. Thus, the ensemble-averaged color charge 
densities vanish, 

{J^A^Mxh'^^Mix)} = , (2) 
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where the ■0ai are the quark fields and A"j, the color generators. As a consequence, 

nR^nc^riB , (3) 

with Ha = {T.^'4'ai{x)^°'4'aiix)) , and 

Uab = {'S,ill>ai{x)-f"i>biix)) = (4) 

for a ^ b. Condition along with condition (||), relates /ia to fi as 

f^ = ^IJ'a ; (5) 



equivalently we write, 



/ia - f • (6) 



Then 

= (7) 

a 

For a 7^ 6, we introduce the chemical potentials /iaf, associated with Hat in such a way that the grand-canonical 
Hamiltonian reads K = H — Yl,ab^^a-btJ^/'i + Aafc)Ei / '^■^2:?/'fci(a;)7°V'ai(a;)]. Following Iwasaki and Iwado |2|], we can 
regard these jlab as Lagrange multipliers that ensure condition . The hermiticity of K requires flab = i^la ■ 

Cooper pairing between quarks, which we assume to be in a channel with zero total angular momentum, is most 
generally characterized by a ANjNc x ANfNc gap matrix in color, flavor, and Dirac space fl^ , 

A(fc) = A(i)(A:o, + A(2)(fco, |k|)7-k7V + ^^^\k^, |k|)7°7' + AW(A;o, |k|) 

+A(5)(fco, |k|)7.k70 + A(6)(fco, |k|)7.k + A(^)(fco, |k|)7.V + A(«)(fco, |k|)7° , (8) 

where k is the relative four- momentum between the two quarks forming the Cooper pair, k = k/|k|, and A^"^]^ are 
the NfNc X NfNc matrices denoting the pairing of the quark of color a and flavor i with that of color b and flavor j. 
The gap A is formally related to the pair amplitude {^p^ {x)'ip{y)) via p^j2^ , p^ 

f J4^^ „.\„tQ(x-y)^ai^°' f 



n odd J y ^) J 



where ip^^^ = CtJ)'^^ is the charge-conjugate spinor (C = 17^7" in the Pauli-Dirac representation), F^^^^ is the full 
quark-quark-gluon vertex, F*^^^^ is the full antiquark-quark-gluon vertex, D is the full giuon propagator, and the 
Matsubara frequencies are given by = imrT. (In Sec. Ill we analyze this gap equation in the weak coupling limit; 
see Eq. (|53|).) The Pauli principle requires that A^"' satisfy 

A%{ko,\k\) = Ai%{-ko,\k\) , n=l,...,6, (10) 



Al^(A:o, |k|) = -Ai:>,{~ko, |k|) , n = 7,8. (11) 

For massless quarks the condensates are eigenstates of chirality |^^. Wilson renormalization-group analyses jlj] 
show that the pairing instability occurs between quarks of the same chirality, rather than between the left- and right- 
handed quarks. We thus ignore the terms A*^^^, A^^^ A*^^^, and A*^*^ in Eq. (^), which are associated with quarks of 
opposite chirality |9|. The A^"*) and A^^) terms in Eq. (^ correspond to condensation in the odd-parity channel . 
Effects of instantons, which prefer even-parity to odd-parity condensates ||l^, lead us to drop these terms. Since only 
A^^' and A*^^) remain in Eq. (||), it is convenient to introduce 

0±(A;o, |k|) s AW(fco, |k|) T AW(fco, |k|) , (12) 
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and rewrite A(fc) as 



A(fc) = ^'[^+{k^, |k|)A+(k) + (/)_(fco, |k|)A-(k)] . (13) 

Here, the 

A±(k) = p- (14) 

are energy projection operators for nonintcracting massless quarks, and (/)± denotes the quark-quark (antiquark- 
antiquark) pairing gap. 

We divide the difference Ar2[0±(A;o, |k|);/iah] of the thermodynamic potential near Tc of the superfiuid phase, ilsF, 
from that of the normal phase, f^N, into two parts: 

= Oo + f^cN , (15) 

where 



17o = Ar![0±(fco,|k|);/i„b = O] (16) 
has the form of a usual Ginzburg-Landau free energy, and 

17cN- Ar![0±(fco,|k|);/i,b]-f7o (17) 

is the correction to fio induced by the color chemical potential differences in the superfiuid (the normal phase satisfies 
jlab — because of rotational invariance in color space). In expressing J7o in terms of we follow the line of argu- 
ment of Mermin and Stare |l^. Note that the grand-canonical Hamiltonian, described by the chemical potential, /i/3, 
common to all combinations of spins, flavors, and colors, and zero masses, keeps ilo invariant under global U(l) gauge 
transformations and rotations in color and flavor space; under global f7(l), and special unitary color and flavor rota- 
tions of the held operators, e^'^'UcUftp, the (/)± transform as {(l>±)abij e' mj- 
Taking into account the condition ( p^ imposed by the Pauli principle, we obtain, up to 0(A*), 

+a-TT{q}U^)F + /3r[Tr(0l(/._)i.]2 + l3^Tii{<j)U_)^]F . (18) 

Due to the relation, A^A^ — 0, the contributions of the quark-quark condensates separate from those of the antiquark- 
antiquark condensates. The subscript "F" denotes the pairing gap evaluated for the quark or antiquark quasiparticle 
momenta |k| equal to the quark Fermi momentum, fc^?. 

In general, the energies of antiquark excitations, ~ I (|k| + /c^)^ -I- 4>_4>- ) i are much larger than the energies of 

quark excitations, ^ I (|k| — fcj?)^ -I- 4>\(j)+ ) . Thus in the temperature region near Tc where /i/3 is large compared 

with the magnitude of the energy gap, the amplitudes for virtual excitations of antiquark quasiparticles with |k| near 
kF are smaller than those for excitations of quark quasiparticles by a factor ^ O(|0-|_(T Tc)\/ as a consequence, 
and [3^ dominate over a~ and (3'^ , and we set a~ — [3^ = 0. 
The next step is to express the correction term i^cN in terms of and jlab- We need retain contributions up to 
second order in jlab- Since — --dflgp/dfia and Uab — —9flsp/dp,ba for a ^ b, the conditions for color neutrality, 
Eqs. (0) and (0), can be written as 



d/\n dAfl dAfl 



djj.RR djlcG djjiBB 

and 

9Af2 „ 



djlab 



(19) 



(20) 



for a^h. 

Up to 0{ilab4'\) and 0(//^j), the only terms in ficN that are invariant under global U(l) gauge transformations and 
special unitary transformations in flavor space are proportional to TvjCi"^, Ti{(f)^_^(j)j^M) f ^ and Tr(9i)+(/)^A^)i?, where 
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■^abij = Sijflab (21) 



in deducing this structure we have used the relations Tr{(j)ji-(j)\^M'^)F = Tr((/)^0+yU)i7 and Tr{(j)^_^(j)^M'^ 



)F = 



Tr((/)+0^A^)F, derived from {(j)'^)F — (4>+) f (se e Eq. (|10[)). The term proportional to Tr(0+(/)^_A^)i?, which does not 



occur in the weak coupling limit (see Eq. (121)), always has a zero coefficient, even in the strong coupling regime. 
The reason is that this term, of second order in the gap, is not affected by the dependence of the pairing interaction 
on the pairing gap, which can only induce terms of at least fourth order in the gap. We thus write 

= crTrTW^ ^ ^Tj:{(I)14>+M)f ■ (22) 

The first term on the right side of Eq. ( p^ ) comes from the change in the normal-fluid free-energy for non-zero jOi, 
while the second is the correction to the £'((/'+) term in JIq. Because of color neutrality, a and x appear, as we shall 
see, only in the dimcnsionlcss combination in the energy gap and the condensation energy. 

In the remainder of this section, we consider condensates that are antisymmetric in color and flavor. In these 
condensates. Cooper pairing is characterized by the products between the color-antitriplet states [{\RG) — \GR))/V^, 
{\GB) - \BG))/V2, {\BR) - \RB))/V2] and the flavor-antitriplet states [{\ud) - \du))/y/2, {\ds) - |sd))/\/2, {\su) - 

\us))/V%- _ . rn. . 

For total angular momentum J = 0, the Pauli principle constraint (nO) requires that antisymmetry of pairing in 
color space be accompanied by antisymmetry of pairing in flavor space. Such pairing states can occur in the weak 
coupling limit, because one-gluon exchange in the color antitriplet channel is attractive j^]. The gap is characterized 
by the ansatz, 

{4>+)abij = ^ijl^abcAlc , (23) 

where A = (dj^, dg, dg) is a general complex matrix. Aia^ — {Aa)i^ denotes the gap for pairing between quarks of 
colors h and c with a ^ h ^ c and of flavors j and I with i ^ j ^ I. 

Substituting Eq. ( [2^ ) into Eqs. (|l^) and (|2^), we derive the Ginzburg-Landau form of fio and ficN in terms of A 
or da: 

= aTT{A^A)F + /3i[Tr(AtA)f ]2 + f3^TT[{A^ Af]F 

= aA + (/5i+/32T)A^ (24) 

with 

a = 4a+ , (3i = 16/3+ -f 213+ , /32 = 2(3+ , (25) 



A^^|d„||, T^^^|d:.d,||, (26) 

a ab 

and 

!^CN = 3a^|/iafc|'-2x^(d:-db)F/iafc . (27) 

ab ab 

One can readily show that under a global U{1), color, and flavor rotation of the flcld operators, -0 — *■ e'^'^UcUftp, A 
transforms as Aia — ^ 6 '^^'^ {U c) abiU f )ij Ajiy . The second order and fourth order terms in A^ included in flo; Eq. ( |24| ) , 
are the only invariants under these transformations. In Ocn, Eq. (p7|), the term of linear order in flab is not affected 
by global U(l) gauge transformations, da e~^"^da, and flavor rotations, da — > Ufda- 

Note that T, which is dimensionless, ranges from 1/3 to 1. For T = 1/3, the configurations for da are determined 

by 

d^.dG = d^.dB = d^.dfl = 0, |dfl|2 = |dG|' = |dB|' ; (28) 
for T = 1, the vectors d/j, dc, and d^ arc all parallel: 

dfl II da II Ab . (29) 



In the weak coupling limit, which we calculate in Sec. IV, the coefficients a, /9i, (32, cr, and x reduce to 
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47VW3)ln( J) , (30) 



/3i = J^^W3) , (31) 



8(^r,)2 

7C(3) 



N{p/3) , (32) 



" 8(^r,)2 

a = -Ni^I/3) , (33) 
X--lnf— )iV(M/3) , (34) 

where 

is the ideal gas density of states at the Fermi surface, and the zeta function C(3) — 1.2020 . . .. 

Two effects beyond weak couphng must in general be taken into account. The first is the modification of the pairing 
interaction due to the pairing gap, which modifies the coefficients /3i and (32 of the fourth order terms in A. The second 
is radiative corrections by the normal medium pl| ] - i.e., quark self-energy, gluon polarization, and quark-quark-gluon 
and three-gluon vertex corrections - which modify the coefficients Eqs. (p0|)-(|34|), mainly through their dependence 
on N{fi/3) and T^. 

We proceed to minimize the Ginzburg-Landau free energy, Eqs. (|4|) and (|27|), with respect to the {da)F- To 
elucidate effects of the color neutrality, it is instructive to start with the optimal expressions for the energy gap and 
the condensation energy in the case in which the flab Sind hence A4 and f^cN vanish. From fig, written up to second 
order in A, we find that thermodynamic stability requires P = Pi + [32^ > for 1/3 < T < 1. It is straightforward 
to show that in this stable region, only two phases occur: T = 1/3 for /32 > and T = 1 for j32 < 0. We discuss the 
physics of the region where (3 < below. 

The order parameters in the T = 1 phase satisfy condition (|2^). All such order parameters lead to states degen- 
erate in energy. Note that in this order-parameter set, any state is identical to an isoscalar, color-antitiplet state 
characterized by {da)i oc Sis (s, the strange flavor), or equivalently, 

{4'+)abi] ^ £abceijs{dc)s , (36) 

to within a constant phase factor and a special unitary transformation in flavor space. In the corresponding condensate, 
the u and d quarks are paired in an isosinglet state. We shall refer to the T = 1 phase as the isoscalar phase, even 
though it contains order parameters having the other orientations in complex flavor space than the s direction. The 
magnitude of the gap in the isoscalar phase is, 

with condensation energy, 

An = ^— ^ . (38) 

The order parameters in the T = 1/3 phase fulfill condition (p8|), and lead to degenerate states. This phase is 
the color-flavor locked phase [ p^ ; the condensate in this phase is characterized by its symmetry under simultaneous 
exchange of color and flavor. States belonging to this phase transform into one another under global U(l) gauge 
transformations and special unitary transformations in flavor space. The simplest among these states is described by 
{da)i OC Sai and (dfl)„ — {dG)d = {dB)s = ka; the corresponding gap matrix is given explicitly by 



'■'+)a 



)abij = HA{5ai5bj - SajSbi) ■ (39) 
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For the color-flavor locked phase, we obtain the magnitude of the gap, 



and the condensation energy. 



a' 



In the original definition of the color-fiavor locked state, in addition to the gap, ka, in the color and flavor 
antisymmetric channel, a gap matrix Ksi^aiSbj + SajSbi) arises in the color and flavor symmetric channel. In the weak 
coupling limit, the gap ka can be generated by the attractive one-gluon exchange interaction in the color- antitriplet 
channel. Nonzero values of K5, on the other hand, are not driven by the one-gluon exchange interaction since it is 
repulsive in the color symmetric channel; they are 0{g) correction to ka that ensures the existence of the ka ^ 
solution to the relevant weak coupling gap equation (see Ref. This implies that the gap ks can be ignored in 

the temperature region near the onset of the pairing. 

We now construct, in the P1-P2 plane, the phase diagram exhibiting the more stable pairing state of the isoscalar, 
and color-flavor locked phases. We first address the question of whether or not Tc, the temperature at which the 
pairing instability of the normal phase occurs, is the same for these phases. The onset of this instability is controlled 
solely by the pairing interactions in the normal phase between quark quasiparticles with zero total momentum. The 
corresponding amplitudes depend on the color, flavor, and quantum numbers such as the total angular momentum, 
chirality, and parity of the quarks involved in the pairing j^, and for the two phases considered here, the quark-quark 
pairs have antisymmetric structures in color and flavor, the same chirality, even parity, and J = at the onset of 
the pairing instability. We thus conclude that no difference in Tc arises between these phases. The amplitudes in the 
normal system driving the instability do not distinguish between the final possible paired states. 

Figure 1 shows the map of the isoscalar and color-fiavor locked phases in the /3i-/32 plane, for the case in which Ai, 
or equivalently, F, is set equal to zero. The regions subtended by these phases are restricted by the region in which 
/3 > 0, for 1/3 < T < 1, is violated. By comparing the condensation energies Eqs. ( |38| ) and (p]), we find that in 
the weak coupling limit, described here by Eqs. (|3^)-(|3^), the color-fiavor locked state is favored over the isoscalar 
channel, a result consistent with the conclusion drawn from the weak coupling analyses |0,[l4] at zero temperature. 

















CFL 









/ 

IS 


First order 
transition 









^1 



Fig. 1. Phase diagram in the Ginzburg-Landau regime, showing regions where the isoscalar (IS) and color-fiavor 
locked (CFL) phases are favored, when the constraint of color neutrality does not affect the free energy, F = 0. The 
(3i, are the fourth order coefficients in the Ginzburg-Landau free energy, Eq. (p4|), and F (Eqs. ( ^ ) and (|27|)) 
describes effects of the constraint of color neutrality. The cross indicates the weak coupling limit. In the region of 
first order transitions, the overall fourth order coefficient, /3, is not positive definite. 
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The region in which /3 < corresponds to a first order phase transition from the normal to the superfluid state 
at a temperature greater than at which a{T^) vanishes, and the normal state becomes unstable against Cooper 
pairing. The situation is similar to condensed matter systems in which the Ginzburg-Landau free energy contains 
a term cubic in the magnitude of the order parameter, and also similar to the situation in which the chiral phase 
transition in QCD with two-flavor massless quarks changes from second order to first order at a tricritical point in 
the iJ, versus T plane To see this structure, we add to the Ginzburg-Landau expansion up to fourth order in A, 
Eq. (p^), the sum of the sixth order invariants, 

= (7l+72T-t-73e)A^ (42) 

where A and T are given by Eq. (p6|), and 



db)F{dl-d,)F{d:-da)F . (43) 



Let us assume that the overall coefficient of , 71 -f 72T + 73^ = 7, is positive. Then, the local minimum of fig + is 
reached at a nonzero value of A that can be calculated as Aq = |/3|/37+ [(/3/37)^ — a/37]^/^. The critical temperature 
Tc can be determined from the condition r2o(Ao) -I- 06(Ao) = 0, or equivalently, a(Tc) — /3^/47, where we can ignore 
the temperature dependence of (3 and 7. Note that T^. is greater than T°. We thus find that Aq = |/5|/27 at T = T^, 
i.e., the pairing gap is discontinuous at the transition point. We remark that this argument is only applicable to the 
case in which Tc — ^ . In the case that /3 < and 7 < 0, one must go to higher order to determine the critical 
temperature and the discontinuity of the order parameter at the transition. 

We turn to the Ginzburg-Landau energy, Eqs. ( p^ and (^), for a color-singlet system. As we shall see, the 
constraint of color neutrality acts to modify the pairing gap without removing the degeneracy of the order parameters 
occurring when the differences in color chemical potentials vanish. Generally, the deviations of the chemical potential 
differences flab from zero can be determined from the color neutrality conditions (|l^) and ( pO[ ) as 

nab^^ma-dl)F-5abY.\dc\l]. (44) 

c 

Substitution of this expression into Eqs. (^ ) and (|2^) leads to the Ginzburg-Landau form, 

Ar! = aA+[(/3i-r) + (/32 + 3r)T]A2 , (45) 

with 

'^^ . (46) 



cr V3 

The coefficient F is positive definite as long as the color-singlet system is thermodynamically stable against color 
fluctuations, i.e., {d'^An/dilabdij-ba){d^)F = 6(t < 0. This condition holds in the weak couphng limit, Eq. (|3^); 
vacuum polarization effects in the low-density regime, which lead to color antiscreening, or positive color susceptibility, 
—d^ft/dfLabdflba > 0, are in the direction to preserve this condition. We thus assume F > for arbitrary coupling 
constant g. 

Note that Eq. ( ^ ) is identical with Eq. ( p^ ) except that (3i and P2 are replaced by /3i — F and /32 + 3F, i.e., inclusion 
of the color neutrality results in a renormalization of /3i and /?2. This is due to the fact that the constraint of color 
neutrality itself is invariant under color rotations and thus chemical potential differences between colors yield the 
fourth order terms invariant with respect to rotations in color space. 

We then obtain from Eq. ( ^5|) the two sets of optimal order parameters, T = 1/3 (the color-flavor locked phase) 
for /32 + 3F > and T = 1 (the isoscalar phase) for /32 + 3F < 0. The resultant magnitude of the gap in the isoscalar 
phase can be written as 

^ " "2(/3i+/32 + 2F) ' ^^^^ 



with condensation energy. 



An = . (48) 

4(/3i + /32 + 2F) ^ ^ 
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We observe from the energy gap and condensation energy, Eqs. ( [47D and that A and |Afi| are suppressed at 

fixed Pi and P2 by the parameter F characterizing the color neutrahty. This suppression, together with 7^ as can 
be found from substitution of condition d29) into Eg. (Ej), comes from the fact that the condensate is intrinsicaUy 



anisotropic in color space, as we shall see in Sec. [II A. We remark in passing that in the weak coupling limit, where 
only contributions of leading order in ZTd jj. remain, T is dominated by /3i and [32 and hence does not significantly 
affect the gap. In the color-flavor locked phase, on the other hand, the gap size and the condensation energy of the 
color-flavor locked phase are still given by Eqs. ( ^ ) and (p]). For this phase, we can find = from substitution of 
condition (^ ) into Eq. (^) . This vanishing reflects the fact that the condensate in the color-flavor locking is isotropic 
in color space, as suggested by condition (pq). 

We conclude this section by examining the influence of the requirement of color neutrality on the phase diagram 
in the /3i-/32 plane. As we see in Fig. 2, the phase diagram for F > 0, this requirement shifts the isoscalar phase by 
F in the /3i direction and — 3F in the [32 direction. With increasing F, this phase, in which the order parameters are 
anisotropic in color space, moves farther away from the weak coupling point, (/3i,/32) given by Eqs. ( pi] ) and (|3^). 
With increasing F, the color-flavor locked phase begins to cover the region occupied by the isoscalar phase. 

To answer the question of which phase is more stable at low densities will require understanding in detail how strong 
coupling effects, which develop as the baryon chemical potential /i is lowered, modify (32, and F with increasing g. 
These effects are too uncertain for us to be able to predict in general the more stable state near Tc- 




r>o 



Fig. 2. Same as Fig. 1 with the constraint of color neutrality, described by F > 0. 



III. GAP EQUATION 

We turn now to deriving the high-density asymptotic form of the thermodynamic potential difference Af? near 
Tc, as we used in the previous section. The gap equation in the weak coupling limit provides information on the 
overall temperature dependence of the gap and the quasiparticle structures, in contrast to the Ginzburg-Landau 
approach developed in the previous section, which concentrates on limited temperature region near Tc and coarse- 
grained features of the superfluid. For the moment, we consider an arbitrary color neutral condensate of the form 
(|13|). For the two optimal pairing states considered in Sec. ||, we obtain the behavior of the energy gap from the 



finite-temperature gap equation, with the constraint of color neutrality. In Sec. IV, we derive the coefficients given 
by Eqs. (pO|)-(^^ from the structure of the energy gap. 

In obtaining the gap equation in the weak coupling limit, it is convenient to first introduce the notation for the 
quark field, (V'ai, V'S)- We then write a self-consistent Schwinger-Dyson equation to determine the proper self energy 
E(/c) up to Cl(ff2) 
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n odd 



7''A"/2 



-7^(A")^/2; -ri^^Vf^^ 

where the summation is over Matsubara frequencies, go — imrT, with n odd for fermions, 

= l^G(2i)(fc) G(22)(fc) j - 1^ A(/fc) 7fc-70X^ 







(49) 



(50) 



is the quark propagator with normal state Hartree-Fock contributions ignored, A = 7'^A^7'^ and M.abij ~ ^ij{^abi^/i~ 
p.ab): and D'^j^ is the gluon propagator, specified below. 
The self-energy and quark propagator are related by 



where 



r (k\- f ^^ + ^"-^ 

is the noninteracting quark propagator. 

The off-diagonal component of Eq. (Bof) yields the gap equation in the weak coupling limit, 



(51) 



(52) 



•' ^ n odd 



2 



where 



G(21) (q) ^ (75 - 7"X^)-1 A(q) \A{q)i^q - j^M^r'^ ~ (79 + M) 



(53) 



(54) 



is the 21-component of G{q). The summation over n on the right side of Eq. ( p3|) leads to the self-consistent gap 
equation. 



A(a-±(k),k) = -lg2^,7/ ^i?^.^(£±(k)-r;£,(q),k-q) 



X7' 



,7=± 



(27r)3 

A(w^(q),q)e,7^(q) tanh 



^.)(q) 

2r 



where the (on-shell) gap matrix A(fc) is given by, 

A(yfc) = A-(k)7°A(fc)A+(k)7° , 

with the Hermitean frequency matrix 

w±(k) = ±e±(k) - ilmM , 

and 



(55) 

(56) 
(57) 



{e.ri{'k))abij = ^((5ac(|k| - ^/3) - Re/iQc)(<^cfc(|k| - fi/3) - Re/Xch) 



(58) 



clfll' 



the energy matrix for quark quasiparticles. Here, we have used A4A4^ = (ReA^)^ + 0{A4'^) and made the approx- 
imation in G^2i)(g)^ Eq. (|4|), that (79 — M'^)^^ A{q){-jq — 7°7W^) w A^(g). This approximation is sufficient to 
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describe the gap matrix up to first order in since ||q| — /i/3||/iab| is much smaher t han t he square of the gap 
magnitude in the momentum region where the gap is appreciable (as we shall see in Sec. [II A ). We have also disre- 
garded the contributions of antiquark quasiparticles to the integral over q in Eq. ( |55| ) ; in the weak coupling regime, 
these contributions are suppressed by one power of A//i compared with the quark quasiparticlc contributions since 
the integral is dominated by the region immediately close to |q| — /i/3. Using the definition (p^, we obtain 

A(fc) = -7^0+(fc)A-(k) , (59) 

and the quasiparticlc energies 

e,(k) = [(|k|-ReA^)2 + 0^(c.,(k),k)0+(c.,(k),k)]i/2 . (60) 
Equation (Bq) can now be rewritten in terms of cj)^ as 



M^±{^),-i^)^-YT.^J (|^^^.^(£±(k)-^e„(q),k-q)Tr[7''A-(q)7''A+(k)] 

xi(A")^0+(^,(q),q)£-i(q)tanh (^^^ . (61) 

The gap equation (^ij) contains the color chemical potential differences flat, which are determined by the color 
neutrality conditions, Eqs. (^ and (^, written in terms of Uat with Ua = riaa] these Uab, characterizing the color 
charge densities, are given by 

"'^'' = ^E/(0ETr[/Gii:)(<z)], (62) 



where 



G^^^\q) = [(79 + I'^M) - A{q){-fq - j°M^)-^A{q)\ ' (63) 



is the 11-component of the propagator G{q). In terms of (f)±, 

G(")(g) = (jq-j^M^) {[{qo + ilmMf - {\q\ + ReMf ~ cl>Uq)cf,^{q)]-^ A+ (q) 

+ [{qo + tlmMf - (|q| - RcMf - cj,liq)<j>+{q)]-' A' (q)} , (64) 

where we have again used the approximation {jq — •j'-' Ai'^)^^ A{q){-yq — j'^Ai'^) « A^{q). As in the derivation of the 
gap equation (^5|), this approximation is sufficient to determine Uat up to first order in Ai. Substituting Eq. ( |6^ ) into 
Eq. (|6^), we find 

nab^TY, f 7^ E 2{{qo-M^-m{qo + tlmMf-{\q\+RcMf~cl)Uq)^-{q)r' 

i ^ ^' n odd 

+ {qo-M^ + m{qo+tlmM)'-i\q\-RcMf-cpl{q)<f>+{q)]-'} . (65) 

J ahii 

The Matsubara frequency summation in this equation yields the usual BCS expression, 

"""^E/ 70{l-<(q)[l-2/(c.+ (q))]+«^(q)[i_2/(-c._(q))] 

i 

+0{4>lct)./^i^)}aU^ , (66) 

where 

is the distribution function for quark quasiparticles, and the BCS coherence factors are given by 



12 



with u^j(q) + w^(q) = 1. In obtaining M up to leading order in 0+, Eq. (|6 



ReA^-|q|]e-i(q) , 
can be simplified as 



(68) 



nab 



(27r)3 



{/(|q|-A4)-/(|q|~A^/3) 



[4(q)/(£+(q)) + ^^+(q)(i - f{e+{m]M=o}aM^ 



(69) 



where we have used the fact that £+(q)|imAi=o — e-(q)|imX=o = e(q) and thus 0+(e(q), q) = (^+(— e(q), q). 

In specifying the gluon propagator D^fjv = fc — g), it is essential to take into account the long-range nature of 
the color magnetic interactions [p^|jl^ , |20| , |34 . Once the effects of a normal medium are included in the random-phase 
approximation (RPA), these interactions are dynamically screened by the Landau damping of the virtual giuons 
involved ^j^. Since the dominant contributions from D{p) are peaked around p = 0, we focus our attention on 
bo I ^ m/3 and |p| ^ ^/3. The gluon propagator in the Landau gauge can be written within the RPA as 



with the dimensionless transverse and longitudinal polarizations, 



m]jXL{pol\v\) ' 



Xt{x) = — 



■In 



(70) 



(71) 



and transverse and longitudinal projection operators, 



1 - - In 
2 



P 



00 



-in. 



, 



(72) 



(73) 



pL 



2 i'A''^ /it^ 



where the Debye mass is given by 



67r2 



(74) 



(75) 



In Eq. (^), the term proportional to P^ describes the transverse or color magnetic sector, while that proportional 
to P^ describes the longitudinal or color electric sector. Landau damping provides an effective infrared cutoff in the 
magnetic sector, ~ (7rm|,|po|/4)^^'^, for |po| ^ |p|- This cutoff dominates the determination of the energy gap; the 
infrared cutoff given by a putative magnetic mass |^o[] makes negligible difference. In the presence of color pairing, 
certain of the magnetic giuons acquire a mass as a consequence of the Meissner effect |3|,Q, while the Debye mass 
of some of the electric giuons is reduced . We ignore here these effects of the superconducting medium, although 
they possibly change the g dependence of the gap matrix in the weak coupling limit. 

For the purpose of solving the gap equation ( |6l| ) to leading logarithmic order in g, it is practical to approximate the 
polarized gluon propagator ([70|); we first set « |pp and xl(x) « 1, neglect the real part of xrix), and replace 
ImxT(a;) = —TTx/4 for small |a;| by —{■Kx/A)9{y/TTmo/'2, — |p|). The effect of this latter replacement is to introduce a 
high momentum cutoff for the Landau-damped magnetic giuons. Then, we replace the gluon propagator by its real 
part, so that 



DfAp) 



-Re 



|p|2 - iTTmlpa9{^mD/2 - |p|)/4|p| 



(76) 



This form is equivalent to that in Ref. |^^. In the gap equation (|6l|), we may replace |q| by |k|, and then replace q 
by k in the explicit Dirac structures Tr[7^A^(q)7''A''"(k)]P^f'^^(fc — q), so that both these terms become —2. We then 
find that the gap is determined by 
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6+(a;±(k),k) 



' 4 



(2^ 



|q-kP 



|q - k|4 + (7rm|,(r;£,(q) - £±(k))/4|q - \^\YB{^inol2 - |q - k|) 
1 



|q — kp + m 
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X T ( A" (t^,, (q) , q)e^ ^ (q) tanh 



^»;(q) 

2r 



A" 



(77) 



where is the angle between k and q. The dominant contribution to the integral in Eq. (77) conies from colinear 
scattering between the quarks of momenta k and q. Concentrating on the corresponding momentum region, cos ~ 1 
and |k| ~ |q|, allows us to solve for the gap to subleading order in g by reducing Eq. ( [77| ) to 



0+(w±(k) 



„2 fca r ^ 



|?7£^(q) -e±(k)| 



xi(A"f 0+(c^,,(q),q)£-i(q) tanh 



Mt 

^>?(q) 

2T 



In 



A" 



(78) 



with Mt = \fTiraDl1. Below, we shall consider the two optimal pairing states analyzed in Sec. ^ and estimate from 
Eqs. (^ and (^ the chemical potential differences between colors and the reduction in the pairing gaps induced by 
imposition of color neutrality. 



A. Isoscalar, color-antitriplet channel 

Cooper pairing in the isoscalar phase is described by a gap of the form (|3^), to within a constant phase factor 
and a special unitary transformation in flavor space. This gap matrix is generally accompanied by nonzero chemical 
potential differences between colors, as we found in the Ginzburg-Landau regime in Sec. ||. We estimate here the 
chemical potential differences, using the gap equation (^sh and the color neutrality conditions (||) and (||) with (|69|). 



Note that the right sides of Eq. (69) as well as Eq. ( |78| ) with ImA^ = are integrals of an even function of e(q), 
since under e — > — e, v^. and / <-> 1 — /. It is thus convenient to diagonalize the matrix £^(q) up to first order in 

the chemical potential differences, entering via M.. We first transform bases in color space from the original (RGB) 
basis to a new basis [R'G'B') to diagonalize the matrix (/)^,_(q)(/)+(q) where (/)+(q) = 0+(£(q),q). When = 0, the 
unitary matrix U that carries out this diagonalization also diagonalizes e^(q)- As we shall see, U is determined solely 
by the direction of the vector d, with components da = (da)s, in complex color space and hence is independent of 
q. In the isoscalar channel where the jj-components of (l)\(j>+ obey ((/)^). (/)+)(„„) = {4)^(1)+) (^dd)^ {4>\4'+){ss) = 0, and 
(0^0+)(jj) — for i 7^ J, it is sufficient to diagonalize the submatrix {4'\<t'+)(uu)- The result reads 

/|d|2 0\ 

u\4,l(ly+)(^,,)u = |d|2 , (79) 
\ 0/ 

where 

M = C/u« = (Ui?',UG/,Us') , (80) 

with 

d*x(e X d) e*xd* , , 

Ufl' = — I , ug' = , ub' = d . (81) 

e X dp V I® ^ '^l^ 

Here, d = d/|d|, and e is an arbitrary complex unit vector that satisfies e x d 7^ 0. Equation ( |79| ) describes the 
pairing state with the convention dj^i = dci = and ds' 7^ 0. The fact that the the matrix {4>+4'+) (uu) is effectively 
of rank two, corresponding to a reduction of the color symmetry from SU(3) to SU(2), implies that the condensate in 
the isoscalar, color-antitriplet channel is anisotropic in color space. 
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The unitary matrix U can be taken to be block diagonal: = U^dd) = f^(ss) = u and U(^ij) = for i 7^ j. Then, 

multiplication of the gap equation (^q) with ImAi = by C/ on the right and C/t on the left gives 



/ 10^ 
wt[^+(k)](„d)M -10 
\ 0, 



^1^1 U 1- 












1 


°) 


Vo 





0/ 



-5- 




I + In 

,0 1 



Xju''{X°')'^uu'^ (q)] uuh^^^^ (q) tanh 




where ej-jj) is the w-component of e, and we have used — S(dd) a-nd (p+^du) = ^4'+{ud)- 
terms can be transformed, via the Fierz identity (A")afc(A")c(j = —{2/3)SabScd + '^^adSbc, to 

^ul^,{X")abUab'U*c'iX")cdUdd' = " ^ (^a'f,' (5c'd' + 2(U6/ -U^' ) (Ua' 'Uc' )* • 

abed 



(82) 

The Gell-Mann matrix 
(83) 



Before taking into account the effects of nonzero , it is instructive to understand the structure of the gap equation 
with Ai = 0. We first choose the complex unit vector e in such a way as to simplify the matrix u^(j)^(^ud)U, and then 
write down a reduced gap equation. Because of the invariance of the grand canonical Hamiltonian under color rotations 
in the absence of A^, it is sufficient to treat the states in which d x d* = 0. These states are related to the states in 
which d X d* ^ by a gauge transformation, d e~'^^'^d, and a special unitary transformation, d Ucd. 

For the states satisfying d x d* = 0, we can write d = |d|d exp(i(^o), where d is the real unit vector proportional to 
d, and (pa is the phase common to all the elements of d. Let us take e to be a real unit vector that satisfies e x d 7^ 0. 
We then obtain 



/ -|d| 0\ 

u'^ <l)+(ud)U = |d| exp(2i(/3o) , 
\ 0/ 

as well as the relation, Ua'-Uf,/ — Sa'b'[Sa'R' + ^a'c exp(— 2i(^o) + Sa'B' exp(2z(/3o) 
reduces to the single equation. 



(84) 



The matrix equation (183) thus 



|d(k)|=-^y^ ^Iqljgln 
x|d(q)|£;"i(q) tanh 



\E\ci)-E\\,)\ 



Af2 



In' 



4|qP 



In 



ml 



4|q|2 



2T 



where 



(|q|-g)' + |d(q) 



1/2 



(85) 



(86) 



1/2 



is the R'R' and G'G' component of the diagonal matrix (^[u^e(„„)(q)w]^^oj • Note that Eq. ( |85| ) with dji = da = 
is equivalent to the gap equation usually analyzed for the isoscalar color-antitriplet channel (see, e.g., Refs. 32|). In 
the weak coupling limit of interest here, where a tiny energy gap arises in the momentum region immediately close to 
the Fermi surface, we can replace the E^{x) (x = k, q) in the logarithmic term of Eq. ( ^5|) with the S'iJ'-component 
(|x| — fJ,/3)^ of the diagonal matrix. We then assume in Eq. ( |8q ) that |q| « fi/3 in the last two logarithms and shift 
the integration variable from |q| to |q| — The resulting equation reads 



|d(k)| 



IStt^ 



rf(|q|-M/3) 



tanhf^Uln 

V 2T y 2 



|q|-M/3)2-(|k|-M/3)2| 



|d(q)| 



(87) 
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with 



256.^ f^^ 



5/2 



the cutoff 5, obeying |d(|q| = /i/3)| <C (5 ^ to_d, is chosen so that |d(|q| > S)\ is vanishingiy smail. Here, we have 
used ml = g^fJ,y6TT^ + Oig^T^). 

Let us now consider a coior-singiet system and ask, in particular, whether or not the degeneracy of the order- 
parameter sets is removed in the weak coupling limit by the chemical potential differences between colors stemming 
from color neutrality (||) and (||). As a first step, we write the solution to the gap equation (p7|), following a line of 
argument of Pisarski and Rischke and express the chemical potential deviations jlab in terms of the obtained gap. 
We then examine how such deviations in turn affect the gap itself. As we shall see, the degeneracy is not removed up 
to subleading order in g. 

We now implement a process developed by Pisarski and Rischke to convert the integral equation ( p7| ) into an 
equivalent differential equation [Eq. ( |92|) below] via the approximation. 



|q|-A./3)2-(|k|-/V3)2| 



0(||q|-/x/3|~||k|-/i/3|)ln 
+0(||k|-M/3|-||q|-/./3|)ln 



With this approximation and the change of variables. 



In 



26^/3 



\k\~f,/3\+E(k) 



, y = In 



26^/3 



tx/3\+E{q) 



, Xi? = In 



q| - /V3| 

bfi/3 
|k|-/x/3| 

2bfi/S 
£^(|k| =m/3) 



we obtain 



|d(.T)| = f dy tanh (^^^ |d(y)| 



dyytanh(^) \d{y)\ 



(89) 



(90) 



(91) 



where g = g/Sy/^ir, and ln(&/i/3||k| — /-*/3|) ~ x and ln(&/i/3||q| — /i/3|) ~ y have been used. Differentiation of Eq. 
( ^ ) with respect to x leads finally to 



''"l'*<"='l.-Snanh('Si),|<iW| 



dx'^ 



2T 



We summarize the results of Eq. 
jap is given by 



), as obtained by Pisarski and Rischke 



(92) 



At T = 0, the magnitude of the 



\d{x)\ = -b^le-^^'Hmigx) 



(93) 



It was pointed out that the BCS-like exponential term and the sinusoidal x dependence arise from nearly static 
magnetic gluons that undergo Landau damping and mediate the long-range part of the magnetic interactions, and 
that both the higher frequency magnetic gluons, which are little affected by Landau damping, and Debye-screened 
electric gluons play a dominant role in determining the pre-exponential factor. 

The overall coefficient of the n/g^ in the prefactor is usually considered to be correct up to a factor of order unity, for 
the reason that radiative corrections such as quasiparticle wave function renormalization and vertex corrections 
p2[ do not modify the asymptotic form of the gap. However, it is still uncertain the extent to which this result is 
exact since contributions of the Meissner effect and Debye screening in the superconducting medium remain to be 
clarified |34). 

Expression ( |9^ ) indicates that the gap takes on a peak value |d(xi?)| = (26/i/3)e~'^/^^ at |k| = /z/3, smoothly 
reaches half its peak value at ||k| — fj,/3\ ~ |d(a;_F)|e'^/'^^, and becomes of order {g In g)\d{x f)\ at ||k| — /L^/3| ~ mn- 
The critical temperature Tc is given by 



Tc = —\d{xF,T = 0)\ 

TT 



(94) 
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where 7 = 0.5772 ... is the Euler constant. It was stressed by Pisarski and Rischke |^ that expression ( |94[ ) is of 
the usual BCS form, a feature stemming from the fact that in the weak couphng hmit |d(a;,r)| ~ \d{x,T = 
0)\{\d{xF,T)\/\d{xF,T = 0)1). This expression, except for a factor of order unity due to quasiparticle wave function 
renormahzation, was reproduced by Brown, Liu, and Ren fs^ j using finite temperature diagrammatic perturbation 
theory in the normal phase. Since A4 vanishes at T = Tc, the critical temperature itself is unchanged by the color 
neutrality constraint. 

We turn to the calculations of the chemical potential deviations ^ab as functions of |d| given by Eq. (|9^); we consider 
only zero temperature for simplicity. Up to lowest order in Ai/fi, the difference between nab given by Eq. ( |69[ ) and 
the corresponding normal-phase value Sabnyi = Sab'^iN f / 3) N (fi/ S) fi/ 3 can be written as 



riab - "N<5a6 = 2N f N {p./3) /lab 



E 

=R',G' 



.,iSn 



(95) 



where 




Sn 



is the number density excess due to pairing |37 
(M) and M) as 



m/3 



E 

q|>M/3 



m/3 



(96) 



and V is the system volume. This excess can be estimated from Eqs. 



37r 



5n^^Nit,/3)\d{xF)\' 



(97) 



Using these expressions in the color neutrality conditions (^) and we obtain, with the help of Eqs. (0) and (^), 

5n 



flab = {Sdadl - Sab) 



97V(/i/3) 



(98) 



The chemical potential deviations flab are thus of order \d{xF)\'^ /g^j,. 

Up to subleading order in g, such chemical potential deviations do not modify any of the factors in Eq. (|9^) for the 
magnitude of the gap |d|. This is because ||q| — fj./3\\p.ab\ ^ |d(a;F)p in the momentum region, ||q| — /x/3| < 6, where 
the gap is appreciable. We remark that this robustness of the exponential term and the leading pre-exponential factor 
in Eq. (93) is supported by the renormalization-group analysis developed by Son p8[ |. For example, let us consider 
the case in which — da — and ds 0. In this case, Cooper pairing occurs between the quarks of color R and 
G with momenta ±k close to the Fermi surface having the Fermi momentum kp — /i/3 + flRR — n/S + jlcG- In 
the renormalization-group analysis, the scattering amplitude between these quarks is characterized by the parameter 
t — — In ||k| — kp\. The scattering amplitude becomes singular, as Gorkov originally showed, at the onset of pairing. 
The value of t at which this singularity occurs - the Landau pole ti, - tells us the scale of the energy gap according 
to |d(|k| = kp)\ ~ mne^^'^ . Up to subleading order in g, non-zero fiRR, has no effect on i^,, which behaves as 
tL ^TT/2g + &\iYg + ■ ■ ■. 

In the weak coupling limit, up to subleading order in g^ we thus find that the degenerate set of d's obtained with 
= persists in the color-singlet system studied here, a result consistent with that obtained from the general 
Ginzburg-Landau theory in Sec. ||. Nevertheless, it is instructive to examine the leading contribution of flab to the 
zero-temperature gap, in view of the fact that a strongly coupled quark system may develop a considerable energy 
gap (see, e.g., Ref. |l|]). 

The properties of this contribution can be roughly understood by finding the eigenvalues of £^„„-, (q) , = sfii^) (q) , the 
square of the quasiparticle energy given by Eq. (|60|). This analysis is only applicable to the case in which d x d* = 
and hence ImM — 0. (If ImM 7^ 0, one must go back to the original gap equation ( [78| ) dependent on the matrices 
and LOri, which are not diagonalized by the same unitary transformation.) Up to ©(A^), these eigenvalues are given 
by 





(|q| 


3 


- 6flR 


El,{ci) = 


(|q| 


M 
3 


- 5flG 


EUci) = 


(|q| 


3 


- Sub 
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where 



7r|d(a;F) 



(100) 



is the chemical potential shift from /i/3. The order of the eigenvalues corresponds to that in Eq. (|79[). Note that the 
chemical potential shift 6fi, stemming from the anisotropy of the condensate in color space, is independent of its color 
orientation. This reflects the fact that the properties of a color-singlet system do not depend on the choice of color 
axes. It is also important to note that Sfj, is positive definite. This feature arises due to the property discussed in 
note Not only does positive enlarge the Fermi surface of the gapless quarks of color B', but it also decreases 
the Fermi energy of the quarks of colors R' or G", leading to overall suppression of the energy gap. 

In order to calculate in detail the 0{flab) correction to the magnitude of the gap in the case in which d x d* = 0, 
we note that U diagonalizes A4 and hence e^(q). This is evident from the relation, 



/i 7r|d(a:F) 
3 3gfi 



Sab - 



2gfi 



(101) 



obtained from Eq. (^8|). It is thus straightforward to show that the gap equation ( |82| ) reduces to Eq. 
the quasiparticle energy E{q) is replaced by Ej^r{q), =i?G'(q)- The solution to this equation reads 



in which 



|d(x') 



- exp 



; In 



2Cibfi/3 



\\k\-fi/3 + Sfi/2\+ER,{k) 



where the correction factor is 



Ci 



1-9^, 



(102) 



(103) 



and we have used m| 

Note that Ci is less than unity. As a consequence the anisotropy of the condensate in a color-singlet system acts 
to reduce the gap size. 

This reduction, together with the reduction in the densities of states at the Fermi surfaces for colors R' and G" 
due to S^,R',G'i leads to a decrease in the condensation energy of the system. These features arise despite the fact 
that there are no Fermi momentum differences between the two colors associated with the pairing. This situation is 
in contrast to the case of ferromagnetic superconductors and to the case of color-flavor-locked quark superfluids 
with nonzero strange quark mass ]39| ], where Fermi momentum mismatches between the electron spins and between 
the quark flavors involved in the pairing, respectively, play a role in reducing the magnitude of the gap. 



B. Color-flavor locking 



We proceed to examine color-flavor locking, described by the order parameter (^) to within a constant phase 
factor and a special unitary transformation in flavor space. We first note the relation jlab oc X]i[('/'+<^+)|k|=/^/3]abfj + 
const. X Sab, which can be derived up to first order in from the color-singlct conditions (||) and (^) as well 

as Eq. (|69|). Here, the constant afhxed to Sab ensures Eq. (|^). Substituting Eq. ( |39| ) into this relation we find 
'Yli[{4'^+4>+)\\i\=ti/'i\abii = 4(5Qf,|KAp, and thus obtain Ai = Q. This result, coming from the isotropy of the condensate 
in color space, is consistent with that from the general Ginzburg-Landau analysis in Sec. |l[ We can thus derive the 
behavior of the energy gap at T = and near Tc from the gap equation (|7^) with M. = 0. Substitution of ansatz (|39| ) 
and = into Eq. (fz^) leads to an equation for ka, 



'^A(k) = ^£rf(|q|-Ai/3)iln 



S(8)(q) 



tanh 



;|q|-M/3)2-(|k|-/V3)2| 

^^(8)(q) 



2T 



2KA(q) 
S(i)(q) 



tanh 



g(i)(q) 

2T 



(104) 
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with 



i?(8)(q) = ((|q|-A^/3f + l'^A(q)P 
£^(i)(q) = ((|q|-M/3)' + 4|AJ^(q)| 

Here, we have used the color structure, 



2\l/2 



^ ^ 2 



(105) 



(106) 



and the diagonahzed form of the quasiparticle energy squared, 



(107) 



with 
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{Ru, Gu, Bu, Rd, Gd, Bd, Rs, Gs, Bs) 


in the original 



(108) 



here the bases are taken to be {Ru,Gu, Bu, Rd,Gd, Bd, Rs,Gs, Bs) in the original color-flavor space and 
(R'u', G'u', B'u', R'd', G'd', B'd', R's', G's', B's') in the transformed color-flavor space. 

At T = 0, the gap equation ( |104D is equivalent to that obtained by Schafer ||ll|. The analogy between Eq. (104) 
and Eq. ds^) allows us to write a a in the form (^) valid up to subleading order in g 1ll| : 



KA(a:(8)) = 2 / exp(z(^A)— exp I - — 



where 



^(8) 



= In 



|k|-^/3|+i?(8)(k) 



(109) 



(110) 



and if A is the phase of ka- 

At finite temperatures, the assumption that ^^(k, T) ~ K^(k, T = 0)(Kyi(|k| = fi/3,T)/ ka{\^\ = iJ-/3, T ~ 0)) 
holds as in the isoscalar, color-antitriplet channel. At T = Tc, where £'(8')(q) — E(^i^{q), the gap equation (|l04| ) has 
the same structure as Eq. ( ^7\) obtained for the isoscalar, color-antitriplet channel. Thus, the critical temperature Tc 
equals the result, Eq. (|9j), in the isoscalar channel. Such equality, persisting in the strong coupling regime (see Sec. 
H), stems from the fact that both these types of pairing are induced by the same instability of the normal phase. 



IV. GINZBURG-LANDAU REGION AT HIGH DENSITIES 

We are now in a position to calculate, in the weak coupling limit, the Ginzburg-Landau free energy, (|l5|), the 
difference between the superfluid and normal phases near the transition temperature Tc- As earlier, we ignore the 
normal state Hartree-Fock terms in the Schwinger-Dyson equation ( ^9| ) , and identify the thermodynamic potential in 
the normal phase with that of an ultrarelativistic, noninteracting Fermi gas of Fermi energy /x/3. We thus obtain |2|] 



Ail — Afiideal + ^^loop + i^re 



(111) 



where 



deal 



127r2 



TrA^^ 



(112) 
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is the difference in the ideal-gas contribution between the superfiuid and normal phases; 



n 



loop 



rTr[-G(g)E(g)+lnGo^((Z)G(g)] 



is the contribution, again up to C(y^), of quark loops in the superfiuid vacuum; and 



^rcs — 



T 



E 



odd 



;Tr[G(g)I](g)] 



(113) 



(114) 



resulting from the mean-field approximation adopted in writing Eq. (^), guarantees that the energy gap A 
vanishes at T = Tc- 

We now expand the thermodynamic potential difference (111) with respect to A(q) = A(qo — |q| — /^/3,q) and 
fLab- In calculating fiioop up to ©(A'*) and 0{A'^M), we use the quark propagator G given by Eq. (po|), the proper 
self energy E given by Eq. ( pT| ) , and the noninteracting quark propagator Gq given by Eq. ( |5^ ) . 

The exphcit Dirac structure of G can be obtained from G^"), Eq. (H), G^^i), Eq. iM), 



-1 , 



and 



Thus, 



G(12)(q)= ^A{q){-fq-^"M^)-^A{q)-{jq + -f°M)\ A{q){-fq - -f"M^y^ 
G^'^\q) = \hq - 1°M^) - A{q){^q + ^^Mr'Ail)] • 



(115) 



(116) 



Tr[G(g)E(g)] = 2Tr 



1 



rA(g) 



jq + j^M 79 — 7°A1^ 



A(9) 



1 



A(g)- 



1 



'jq + j'^M '''^^ jq — j^M'^ 



and 



Tr[lnGo-i(g)G(<z)] =Tr 



1 



rA(g)- 



1 



79 + 7OAI 79 — 70yU 



A(g) 



^A(g) 



(117) 



79 + 7°7V1 79 — j^M'^ 



A(<z) 



(118) 



Combining Eqs. (117) and (118) with (113) and performing the Matsubara frequency summation, we finally obtain 

f^loop = <^p + f^{oop- (119) 

where 



(0) 
loop 



d3 



1 

'2 
3 

^8 J {27tY 
- cosh" 



(27r)3 

d3g 



tanh 

tanh 
I - m/3 



m/3 



1 



2T 

|q|-A^/3 



1 



Tr[0t (q)0^(q) 



2T J{\ci\^^Ji|if 



1 



2r ; 2T(|q|-/i/3)2 
xTr[</.^q)</)+(q)0^(q)0+(q)] 



(120) 



is the usual Ginzburg-Landau expansion obtained for = 0, and 



(1) 

loop 



1 

2 

' cosh" 



(27r)3 



/x/3 



tanh 

I - m/3 
2T J 2r(|q| - /x/3) 



1 



2T ){\q\^^i/iY 
1 „.,t 



Tr[0Uq)0^(q)A^] 



(121) 



20 



is the 0{A4) correction to the first term on the right side of Eq. ( |12C| ). We have again ignored the contributions of 
antiquark quasiparticles, and used the relation Tr(0^0+7M) = Tr(</i+(/)^7W"^), coming from (j)^{c[) — 0+(q) [see Eq. 
®]- 

We turn now to calculating o|°^p and ^[Hp for the isoscalar, color-antitriplet channel considered in Sees. || and III, 
and derive the high-density Ginzburg-Landau form of Afl used in Sec. ||. We remark in passing that the calculations 

of J^loop for the color-flavor locking reproduce the values of a, /3i, and P2 that will be derived below for the isoscalar 
channel. 

To derive the coefficients of the Ginzburg-Landau free energy for isoscalar, color-antitriplet pairing, we first carry out 
the integrals over q in i^lo^p, Eq. (120), and Eq. (121). Here we evaluate (/)+(q, T), whose q and T dependences 

are effectively decoupled for T <Tc (see note Q), as 

0+(q,T) ~ 0+(|q| =A^/3,r)sin(5y)|T=o , (122) 
where y is given by Eq. (|90|). We also use the approximation, 



tanh 



2T 



{n\d{xF,T = 0)\ - ||q| -^i/3|)tanh 



\q\-^l/3\-K\dixF,T^0)\)- 



2T 
m/3 



m/3| ' 



(123) 



where k is a positive number of order g ^, and xp is given by Eq. (|9^). With this choice, tanh((|q| — /i/3)/2T) ~ 
(|q| - M/3)/||q| - m/3| for ||q| - > K\d{xF,T ^ 0)|, and sin(gy)|T=o ^ 1 for ||q| - m/3| < K\d{xF,T ^ 0)|. We 
thus obtain 



(0) 
loop 



-iV(/i/3) In 



e-y(26/x/3)i/'|d(a;f,r^0)|i/2 



Tr( 



(124) 



and 



^^loop = - In ( — 1 N{p/3)TTicj,l^+M)F ^ xT<^\^+M)f 



(125) 



where the coefficients of Tr(- • ■)f include the leading contributions with respect to g, and the temperature has been 
set equal to the critical temperature Eq. (p3), except in the coefficient of Tr((/)^(/)+)i?. The term fires, expressed as 



fir, 



N{p/3) In 

7C(3) 
8(7rT,)^ 



e^(26/i/3)i/2|d(a;F,T = 0)|i/ 



21 



Tr( 



N{^i/3)TY{cj)l<j>+cl)l<j>+)i 



(126) 



acts as a counterterm to the coefficient of Tr((/)J^0_|.)i?, leading to Ail = at T = Tc 
fires reduces to the usual Ginzburg-Landau part f2o. 



Then the sum of fi 



(0) 
loop 



and 



flo = Af(M/3)ln(^)Tr(0;^ 



-)f + i^^||^iV(A./3)Tr(0^^+^;0+) 



a+Tr( 



^+(P+)F 



P+TT[i<j>lq^+y]F . 



(127) 



From Eq. (|125[) we derive the value of x> (p^)- Equation (127) indicates that — 0, and implies the coefficients in 
Eqs. @)-^ 

The coefficients in Eq. (127) agree with those obtained by Bailin and Love Q for a BCS-type short-range pairing 
interaction and for Ai = 0. The agreement with the term of fourth order in 0+ arises because the main contribution 
to the corresponding integral in Eq. (120) comes from the momentum region ||q| — fi/3\ ^ K|d(a;i?,r = 0)|, where 
the gap is almost flat. The agreement for the second order term in 0+ is obvious. On the other hand, the coefficient 



21 



of Tv{(l)l(j)+M)F in as can be seen from Eq. ( p^ ), is different from -{6/ fi)N{fi/ 3) ln{ A/ T^) with the ultra- 

violet cutoff A as in a BCS superconductor, due to the behavior 0+(q, T) cx sm{gy)\T=o induced by the long-range 

dynamically screened magnetic interactions. 

The difference in the ideal-gas free energy of the superfluid and normal phases, ASljdeai, Eq. ( |ll2|) , gives rise to a 
term proportional to TrAi^, with coefficient given by Eq. (|33|). The color neutrality conditions, (|19|) and (|20|), thus 
imply 

1 /ST \ 

Ma6 = -:^ln — [3da{xF,T)dl{xF,T) - \d{xF,T)\^Sat] , (128) 



3/i \ /i 

in agreement with Eq. (E3) with {da)i — 5isda- Comparing this expression to the T = result Eq. (S^), we find that 



/2ab/|d(a;i?, T)\ is identical at both T = and T ~ Tc- Expression ( 128 ) can also be derived from the same analysis 



that yields Eq. (Bq) if we replace, in Eq. (Dq), 5n as given by Eq. (|96D with 



(129) 



The equality in iiab/\d{xp, r)p in the T = and T ~Tc cases arises from the fact that both have nearly the same 
momentum dependence of the gap. 

V. CONCLUSION 

In this paper we have laid out the Ginzburg-Landau structure of superconducting quark matter. Even in the 
homogeneous case considered here, many questions remain. First, electrical charge neutrality, in addition to color 
neutrality, should be duly taken into account. The system considered here is composed of three flavor massless quarks 
in flavor equilibrium, as characterized by Eq. ([^). Such matter, when normal, is electrically neutral in itself. Below 
Tc, however, it has nonzero net electric charge, unless the order parameter is isotropic in flavor space. This isotropy is 
retained by the color-flavor locked phase; however, in the isoscalar channel, strange quarks remain gapless, leading to 
a deflcit of the total number of strange quarks relative to the total baryon number, and hence to positive net charge. 
This charge, in neutron star matter, would be neutralized by charged leptons. 

It is straightforward to extend the formalism obtained in the absence of leptons to the situation of electrical charge 
neutrality. When the system is in overall beta equilibrium, the quark chemical potentials, for each color a, obey 
^■au + Me = Mad = Mas, whcrc Me is the elcctrou chemical potential. Then the quark chemical potential, /ia, for color 
a, defined by 

Ma = Mai + ftMe , (130) 

is the same for all flavors i; here, qi is the electric charge of the quark of flavor i. This relation replaces condition (|^). 
As a consequence of electrical charge neutrality, Eq. (P) remains valid. We keep the same definition (H) of the fiaa- 
The matrix M., given by Eq. (pT]), can now be rewritten as a traceless matrix 

Mabij = %(Mafc - qiiJ^e^ab) ■ (131) 



Near Tc, Af2 is still given by Eqs. (|2J) and (|2J), since the lepton pressure, which is of order /ig [~ C'(MQb)]' negligible 
compared with I^cn which is of order M^Mah- (Up to leading order in d^, we obtain /ig = — (x/2(t) gi[(d*)i(da)i]F- 



Oiii^l li^) from the electric neutrality condition, dAft/d^ie = 0, and Eq. (44) from the color neutrality conditions 
( |l9| ) and (^). Thus, Me ^ 0{jlab)-) Since electricity does not distinguish between colors, the presence of Me in the 
isoscalar channel [(da)^ = Sisda] leads to a Fermi momentum mismatch between u and d quarks, and thus results in 
a uniform suppression in the magnitudes of the gap and condensation energy; in Eqs. (^ ) and (|4^), /3i is replaced 
by (3i + r/2. Accordingly, in the phase diagram illustrated in Fig. 2, the isoscalar phase is shifted by r/4 in the /3i 
direction and — 3r/4 in the (32 direction. Note that the general flavor-antitriplet states belonging to this phase can no 
longer be reduced to the isoscalar channel by flavor rotations since the value of Me and the gap suppression depend 
on the the electric charges of quarks involved in Cooper pairing. 

A second important problem is how the strange quark mass affects the quark superfluidity. A non-zero strange 
quark mass not only breaks invariance of the grand-canonical Hamiltonian with respect to rotations in flavor space, 
but it also necessitates the presence of a neutralizing gas of leptons even in the normal phase. The resulting 
dependence of the stability of color-flavor locking over the isoscalar state was considered by several authors (see, e.g., 
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Refs. who found that in contrast to the case of the isoscalar channel, color-flavor locking is destabilized by 

the Fermi momentum mismatch between s quarks and u, d quarks produced by a nonzero value of m^. 

It is instructive to consider the effect of in terms of the general Ginzburg-Landau approach constructed here. 
In the highly relativistic regime (m^ <^ /^/3), in addition to the terms in the massless limit, AO contains a term 
proportional to Tr{(j>\_(j)+M^)p, where Mabij — SabSijUii with m„ = nid = 0. This term tends to suppress the pairing 

gap just like the term proportional to Tr((/)^(^_|_A^)F coming from color neutrality. Note that because of color-flavor 
locking, the resulting anisotropy of the condensate in flavor space acts to fix its color orientation. The transition 
temperature Tc is also reduced in the color-flavor locked state, leading to a Tc in the color-flavor locked state smaller 
than Tc in the isoscalar channel, which is independent of m^. The equilibrium phase diagram for the superfluid 
transition, in the /i versus T plane, is sensitive not only to the effects of the color neutrality and the strong coupling, 
as stressed in the present work, but also to effects of electric neutrality and finite strange quark mass. 
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